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Abstract 

Recently reported tentative evidence for a gamma-ray line in the Fermi-LAT data 
is of great potential interest for identifying the nature of dark matter. We com- 

t""""- pare the implications for decaying and annihilating dark matter taking the con- 

straints from continuum gamma-rays, antiproton flux and morphology of the ex- 

CN cess into account. We find that higgsino and wino dark matter are excluded, also 

for nonthermal production. Generically, the continuum gamma-ray flux severely 

> 

• »-j constrains annihilating dark matter. Consistency of decaying dark matter with 

the spatial distribution of the Fermi-LAT excess would require an enhancement 
of the dark matter density near the Galactic center. 



1 Introduction 



Monochromatic gamma-ray lines have been suggested long ago as signature for pair- 
annihilation of dark matter (DM) particles [I]. During the past decades many DM 
candidates have been discussed, the most popular ones being weakly interacting massive 
particles (WIMPs) in the context of supersymmetric extensions of the Standard Model 
(2]. More recently, also decaying dark matter has been studied in detail [3jj6]. An 
attractive feature of decaying gravitino dark matter is the consistency with thermal 
leptogenesis, contrary to standard WIMP dark matter [I]. 

During the past years the Large Area Telescope (LAT) [7], on board the Fermi 
gamma-ray space telescope, has searched with unprecedented sensitivity for photon 
lines from 30 MeV to 300 GeV. Stringent constraints on decaying and annihilating 
dark matter have been obtained from searches in the energy range 30 — 200 GeV based 
on 11 months of data (8j, and in the range 7 — 200 GeV based on 23 months of data [9j. 
The search region used in these analyses covers the whole sky except for the Galactic 
disk (\b\ > 10°) plus a 20° x 20° region around the Galactic center. A similar analysis 
has been performed independently based on publicly available data corresponding to 



27 months in Ref. 10 . No indications for gamma-ray lines were found. 

In a recent analysis, that is based on optimized search regions around the Galactic 
center, and takes 43 months of Fermi-LAT data into account, a hint for a gamma-ray 
feature in the energy range 120—140 GeV is reported |11|12|. When interpreted in terms 



of a gamma-ray line 12 , the significance obtained from the statistical uncertainties in 
the search regions close to the Galactic center is 4.6o\ The significance is reduced to 
3.3a when correcting for the bias introduced by selecting the search regions. While 



the excess is currently under an active debate 13 -18 , the claim has been further 



strengthened by a recent analysis 19 , which confirms the existence of an excess, as well 
as its spectral shape, with even higher statistical significance than the one claimed in 



Ref. 12 . In addition, indications are found that the excess originates from a relatively 
narrow region of a few degrees around the Galactic center, possibly with a small offset 
within the Galactic plane. Hopefully, the question about the existence of the spectral 
feature and its precise properties will be settled in the near future. 

Assuming that the feature is real, the question about its origin is of great interest. 
While an astrophysical explanation might eventually be identified, a spectral feature 
in this energy range can arise rather generically from the annihilation or decay of dark 
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matter particles. In this paper we compare decaying and annihilating dark matter based 
on two prototype models: decaying gravitinos and wino/higgsino-like WIMPs which 
annihilate predominantly into two W-bosons. In both cases we treat the branching 
ratios into (gravitino) and 77 (wino/higgsino) final states as free parameters to 
account for some model dependence. Note that we do not demand thermal freeze-out 
for WIMPs. Higgsinos and winos can be nonthermally produced in gravitino decays, 



compatible with leptogenesis 20,21 , or, alternatively, in moduli decays 22-24 . 

In the following we shall compare interpretations of the tentative 130 GeV pho- 
ton line in terms of decaying and annihilating dark matter in a sequence of increasing 
assumptions: We first consider constraints from continuum gamma-rays, which are in- 
dependent of charged cosmic rays and the dark matter distribution (Section 2). This is 
followed by a discussion of constraints from antiprotons, which depend on the propa- 
gation model (Section 3). We then analyze the implications of the spatial distribution 
of the Fermi-LAT excess (Section 4) and draw our conclusions (Section 5). 

2 Constraints from continuum gamma-rays 

Since the dark matter particle is required to be electrically neutral, the annihilation 
or decay into photons is typically suppressed compared to channels involving for ex- 
ample electroweak gauge bosons or the Higgs boson, as well as quarks or leptons. The 
subsequent decay and fragmentation of these annihilation or decay products gives rise, 
among others, to an emission of gamma rays with a broad spectrum in the energy range 
~ 1 — 100 GeV below the line. Therefore, one expects that the line flux is accompanied 
by an associated continuum flux with the same spatial distribution, whose strength 
depends on the relative size of the annihilation or decay rates into photons and into 
other Standard Model particles. As compared to the sharp feature resulting from the 
direct decay into photons, the continuum emission has a less characteristic spectral 
shape and is therefore much harder to disentangle from the background gamma-ray 



flux in the Galactic center region 25 . On the other hand, one expects a much larger 
total flux in the continuum compared to the line, and therefore it is important to check 
up to which level a continuum emission is acceptable. Apart from the center of the 
Galaxy, the continuum flux of gamma rays can be constrained, e.g., by the emission 



seen within the Galactic halo 26 , from Galaxy clusters 27 28 , or from dwarf galax- 



ies 29 . In general, a comparison of the various constraints from different targets is 
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affected by the uncertainties related to the dark matter density distribution in the 
different environments. 

We will not enter into the details of this discussion here, but instead follow a very 
conservative approach. Namely, we consider the same data sets for the photon flux that 



have been used in Ref. 12 to search for a gamma-ray line signal. However, instead 
of using only a line as the template, we add a contribution from the continuum flux, 
and test whether the combined spectrum yields a satisfactory fit to the data. This 
strategy has the advantage that it is not affected by uncertainties related to the actual 
dark matter distribution, because both the line and the continuum fluxes are generated 
from the same source. Concretely, we consider a spectral template for the fit given by 

dJ {..„ jp \ , dN EG I - BR, dN? ont \ Q (E\-" 

^ = H ( 7) + ^ + iV 7 BR 7 dE ) +P \k,) ' (2 ' 1} 

where E 1 = rriDM^DM /2), BR 7 = av^/avtot ( T Av) an d iV 7 = 2(1) for dark matter 
annihilation (decay). The fit parameter a determines the strength of the dark matter 
flux, while and 7 model the contribution from backgrounds, assumed to follow a 



power law within the considered energy range (see 13 for possible caveats). For a given 
dark matter profile and search region, a can be easily converted into the dark matter 
annihilation cross section av 71 or partial lifetime t 1u into monochromatic photons. The 
continuum spectrum is in general given by the sum of two components. The first one 
arises from redshifted monochromatic photons with extragalactic origin. For decaying 



dark matter it is given by 30 



dNEG VdmPco E 1 ' 2 

dE y/Tl^iHo/c)^ E: 




(2.2) 



Here Q\ = 1—Qm — 0.73 are the cosmological density parameters, assuming a spatially 
flat ACDM model, VL DM h 2 = 0.112 is the dark matter density, H = /ilOOkm/s/Mpc is 
the Hubble constant with h = 0.704, p c o = 1.05- 10 _5 /i 2 GeV/cm 3 is the critical density, 
c is the speed of light, and the factor = AQ^ 1 J An dQ J los Pd m ( r ) appears because 
of the normalization chosen in (2.1). For the search regions 3 and 4 of fl2| and the 



various density profiles, the prefactor in (2.2) is of the order 0.3 — 0.4. 

The second contribution to the continuum gamma spectrum arises from the weighted 
sum of the spectra produced by annihilation or decays, excluding 77(7//), 

dN 1 f 1 dNl 
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where dNj /dE is the number of photons per energy and per annihilation/decay result- 
ing from each mode, and BRy is the corresponding branching ratio. For the line, we use 
the same shape as in [12] based on the instrument response function of the Fermi-LAT 
detector (FWHM= 0.136), which has somewhat broader tails compared to a Gaussian, 



and also smooth the step function appearing in the extragalactic contribution (2.2) 



accordingly. The continuum flux is determined using PYTHIA 8.1 31 . In order to 
increase the sensitivity to the continuum flux we use the data presented in [12] over 
the full energy range 20 — 300 GeV, taking the more finely binned data where they 
are available (80 — 200 GeV). Lastly, we require that 7 > 2.0 as a conservative lower 



limit for the background power law index 32 , while its normalization is freely varied. 



For performing the fit we use the profile likelihood method as detailed in 12 , with 
95%C.L. statistical errors obtained by varying a (with fixed BR 7 ) and profiling over (5 
and 7 until the TS'-value has decreased by 3.84. We checked that our statistical analysis 



reproduces the best-fit values and confidence intervals in 12 using the smaller energy 
window and a single line up to 10% for regions 2-4, both data samples and all density 
profiles. When including BR 7 in the fit, we use A(TS) = 5.99 to obtain the 95%C.L. 
regions, corresponding to the xl=2 distribution that is expected for two parameters of 
interest^ We will comment on the dependence on the various assumptions later on. 

2.1 Decaying dark matter 

For decaying gravitino dark matter, which we consider as representative example, we 
take the decay channels Zu, W£, hv apart from 71/ into account. In the appendix, the 
partial decay widths are given in the case of bilinear R-parity breaking as functions of 
the gravitino mass, parameters of R-parity breaking, gaugino masses and the higgsino 



mass parameter /i. According to Eq. (5.19), the branching ratio into ju is enhanced for 



hierarchical gaugino masses. For example, in the case of wino NLSP with M 2 — 
one obtains the maximal branching ratio, 

BR max (2.4) 

7 2V2G F m 2 3/2 1 ' 



1 By generating a large sample of data that would be expected in the presence of a true signal, 
we checked that the actual distribution of the profile likelihood ratio is well described by xi=i w hen 
keeping BR 7 fixed, while it is somewhat steeper than Xk=2 wncn including BR 7 in the fit, lying between 
a Xfc=i an d Xk=2 distribution. Therefore the confidence regions assuming xt=2 ma y be regarded as 
conservative. 
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which yields BR° iax ~ 3% for m 3 / 2 — 260 GeV. On the contrary, for higgsino NLSP, 
with fi ~ rri3/2 <C M\, M2, the branching ratio into jv is negligible. 
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Table 1: Partial dark matter lifetime r 7!/ in the monochromatic gamma channel obtained 
from a fit to the Fermi-LAT data in the energy range 20 — 300 GeV in search regions 3 and 4 
of Ref. [12), for the SOURCE (S) and ULTRACLEAN (UC) data samples, for various fixed 
values of the branching ratio BR 7i , = t/t^ u and using the Einasto profile. The columns show 
the line energy that yields the best fit, the lifetimes with 95%C.L. statistical errors, and the 
significance a = VTS (without trial corrections) with respect to the background-only case. 



In order to illustrate the possibilities to explain the tentative gamma-ray line with 
decaying dark matter, we show in Fig. [T] the photon spectrum from gravitino decay, 
fitted to the flux measured by Fermi LAT within the search region 3 of Ref. |12|, for 
three values of the branching ratio into monochromatic photons. The feature in the 
spectrum near 130 GeV can be explained for branching ratios of 10% and 1%. For the 
relatively small value 0.25%, the continuum gamma-ray contribution becomes so strong 
that it would overshoot the measured flux if one would require the monochromatic decay 
mode to account for the feature in the spectrum. In this case, the fit to the observed 
spectrum becomes worse, and requires a somewhat longer total lifetime, such that the 
continuum flux is in accordance with the observed one. The best fit values for the 
partial lifetime, as well as the significance with respect to the background-only case 
obtained from search regions 3 and 4 are shown in Tab. [TJ for various branching ratios. 
The significance of the fit remains nearly constant for branching ratios larger than 1%, 
and then steeply decreases. For a branching ratio of 0.1% or smaller, the dark matter 
component cannot yield a significantly better description than the background only. 

We have also performed a fit where the branching ratio is varied together with 



the parameters a,/3, 7 in Eq. (2.1). The resulting 95%C.L. contours are shown in 
Fig. [2j indicating again that branching ratios below 0.25% are disfavored because of 
the continuum contribution to the photon flux. For comparison we also show the 
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Decaying Dark Matter (SOURCE) 





Figure 1: Sum of monochromatic and continuum photon flux originating from dark matter 
decaying into 71/, Zv, hu, W^i^ with branching ratios BR 7 = 0.1, 0.01, and 0.0025, respec- 
tively. The lower lines show the contribution from dark matter decay, and the upper lines the 
sum of signal and fitted power-law background. The dashed blue line shows the flux when 
taking the gamma-ray line only, plus the redshifted extragalactic contribution, into account. 
The gamma-ray flux measured by Fermi-LAT corresponding to the SOURCE (upper figure) 
and ULTRACLEAN (lower figure) data samples is taken from Ref. [12] . The shaded regions 
correspond to 95%C.L. error bands. 
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Figure 2: Combined fit of the monochromatic line at ~ 130GeV as well as the continuum 
gamma-ray spectrum from dark matter decay. The yellow and orange shaded areas are 
95%C.L. regions for the total dark matter lifetime r and the branching ratio BR(^ — >• ^u) 
obtained from the Fermi-LAT data within search regions 3 and 4 or Ref. |12| , respectively. 
The regions excluded from antiproton constraints at 95 %C.L. based on the PAMELA p/p data 
are shown as blue shaded regions, as well as the expected range for gravitino dark matter 
(red shaded area). Further explanations are given in the text. 
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more stringent constraints resulting from fixing the background parameters ((3,y) = 
{P,y)\BR~ v =i to the best-fit values obtained from the fit using a gamma-ray line only 
(grey dot-dashed line), as well as the constraints obtained when using the smaller 
energy range 80 — 200 GeV (grey dotted line) for region 4. We also note that the 
results are rather insensitive to the precise choice of the relative branching ratios among 
Zv, We, hv. 

Clearly, the limits on the continuum gamma-ray flux obtained in this way are very 
conservative ones, while stronger constraints could be obtained, e.g., by taking photons 
produced via inverse Compton scattering into account (e.g. p\) or extending the energy 
range below 20 GeV. However, this would also make the conclusions more dependent on 
details of the background spectrum and therefore on the astrophysical processes close 
to the Galactic center. In this sense, the limits presented here may be regarded as a 
rather robust consistency check, based on the assumption that the background can be 
modeled by a power law in the energy range 20 — 300 GeV. 

2.2 Annihilating dark matter 

One of the most studied candidates for dark matter is the lightest neutralino within the 
MSSM, for which the annihilation cross sections into 77 and jZ, that yield monochro- 
matic gamma rays, are loop suppressed (331. Consequently, the branching ratios into 
these final states are typically at the permille level, and the continuum flux dominates 
the gamma spectrum. 

In the following we perform first a model independent analysis in order to determine 
the requirements on the branching ratio into monochromatic photons, and then check 
whether some well-motivated scenarios are compatible with these constraints. 

In particular, we first assume that dark matter annihilates into 77 and W + W~ 
final states, giving rise to a line at E 1 = itldm and a continuum photon flux in the 
energy range E < tubm- The resulting photon flux is shown in Fig. |3j for several 
values of the branching ratio BR 7 = av~ n /(av 77 + avww)- Similarly to decaying dark 
matter, we find that the feature in the Fermi spectrum at ~ 130 GeV can be well fitted 
for branching ratios 10% or 1%, while the continuum spectrum becomes dominant for 
0.25%. Note that the photon spectrum exhibits a more pronounced dip between the 
continuum and monochromatic contributions as compared to decaying dark matter, 
due to the absence of the extragalactic contribution. The statistical significance and 
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Table 2: Annihilation cross section in the monochromatic gamma channel obtained from a 
fit to the Fermi-LAT data in the energy range 20 — 300GeV in search regions 3 and 4 of 
Ref. [12), for the SOURCE (S) and ULTRACLEAN (UC) data samples, for various fixed 
values of the branching ratio BR 77 = ov~ n j (cru 77 + avww), and using the Einasto profile. 
The columns show the line energy that yields the best fit, the cross sections with 95%C.L. 
statistical errors, and the significance a = VTS (without trial corrections) with respect to 
the background-only case. The small differences in the statistical significances between dark 
matter decay (see Tab. [TJ and annihilation are mainly due to the extragalactic contribution 
as well as a slightly different continuum spectrum. 

best-fit values of the cross section are given in Tab. [2] for various branching ratios, and 
exhibit a similar dependence than for decaying dark matter. One might notice that 
the significance of the fit always decreases when lowering BR 7 , while for decaying dark 
matter it is slightly better for BR 7 ~ 1%. However, these differences are too small 
to be taken seriously. A common feature is the drastic decrease in the significance for 
branching ratios smaller than 1%. As before, we also performed a fit with variable 
branching ratio. The resulting 95%C.L. regions are shown in Fig. |4| Depending on the 
search region and on the data sample used for the analysis, we find that a branching 
ratio larger than 0.4 — 0.8% is required. For comparison, the constraints obtained when 
using a fixed background contribution (grey dot-dashed line) or when using the smaller 
energy range 80 — 200 GeV (grey dotted lines) are shown as well. The sensitivity to the 
continuum spectrum is essentially lost in the latter case for annihilating dark matter, 
due to the drop of the flux below the line. On the other hand, we checked that removing 
the lower limit 7 > 2.0 on the background slope does not affect the confidence regions. 

As a second generic example, we consider dark matter annihilating into 7Z and 
77, giving rise to two lines at E lt \ = tudm and E 1)2 = tudm — MzK&rriDM)- In 
addition, annihilations into W + W~ and ZZ yield a continuum spectrum. Motivated 
by higgsino dark matter, we fix the relative contribution of the two lines according 
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Figure 3: Sum of monochromatic and continuum photon flux originating from dark matter 
annihilating into 77 and W + W~ with branching ratios BR 7 = 0.1, 0.01, and 0.0025, respec- 
tively. The lower lines show the contribution from dark matter annihilation, and the upper 
lines the sum of signal and fitted power-law background. The dashed blue line shows the flux 
when taking the line only into account. The Fermi data are the same as in Fig. [T| 
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Figure 4: Combined fit of the monochromatic line at ~ 130GeV as well as the continuum 
gamma-ray spectrum from dark matter annihilation into 77 and W + W~ final states. The 
shaded areas correspond to 95%C.L. regions (yellow/orange: Fermi-LAT Reg 3/4 best-fit 
region; blue: PAMELA p/p excluded region) for the total annihilation cross section av and 
the branching ratio av~ n /ov. The prediction from higgsino and wino dark matter (red and 
orange dots) is shown for illustration, as well as the thermal cross-section (straight black line). 
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Table 3: Mass spectrum, annihilation cross sections and denning parameters for two scenarios 
with higgsino (H) and wino-like (W) LSP. Particle masses are given in GeV, and cross sections 
in cm 3 /s. The cross sections are rather insensitive to the other parameters, which we fixed to 
A = 0, Afi 3 = rh = 1 TeV and tan (3 = 7 for illustration. For the wino case the [i parameter 
was chosen such that the neutralino-chargino mass splitting is of the same order than the 
minimal splitting expected from loop corrections |34] in order to yield an estimate for the 
maximal possible gamma cross sections. 

to crt> 7 ^/(2crt> 77 ) = 1.68 as well as ov Z z / www — 0.66. However, we freely vary the 
weighted branching ratio BR 7 = (av- yi + 0.5av lZ )/avt o t- In this case we obtain the best 
fit for dark matter masses in the range 132 — 137 GeV, depending on the search region 
and the data sample, and with similar significance than for one line. The resulting 
95%C.L. regions are shown in Fig. |5j which are very similar to the case of a single line 
discussed before, indicating that the lower limit BR 7 > 0.4—0.8% is rather insensitive to 
the detailed assumptions about the cross sections of the various annihilation channels. 

Finally, we briefly remark that both higgsino and wino-like dark matter can be 
ruled out as an explanation for the tentative gamma line based on the large continuum 
photon flux obtained in both cases. The cross sections for two representative scenarios 
are given in Tab.[3j leading to BR 7 ~ 0.08%(0.2%) for higgsino (wino)-like dark matter. 
Even when allowing for a boost factor or a mixed dark matter scenario in order to 
achieve a line flux of the required size, the continuum flux overshoots the measured 
flux at energies around 20 GeV, which is illustrated in Fig. [6| In particular, we remark 
that this conclusion is independent of the dark matter distribution, and therefore of the 
uncertainties associated with it. According to the recent analysis in (35], electroweak 
corrections could even further reduce the 77 and 7Z cross sections for low-mass wino 
dark matter by a sizeable amount. 
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Figure 5: As Fig.|4j but for a model with two lines from annihilation into 77 and jZ, as well as 
continuum photons from annihilations into W + W~ and ZZ final states. The axes correspond 
to the total cross section av and to the weighted ratio (crt> 77 + 0.5av 1 z)/&v, respectively. For 
the fit we fixed the ratios <tu 7 ^/(2o"v 77 ) = 1.68 and avzz /&vw w = 0.66 motivated by the 
higgsino dark matter scenario. The actual prediction from higgsino and wino dark matter is 
also shown (red and orange dots). 
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Figure 6: Predicted gamma-ray flux for higgsino dark matter (upper graph) and wino 
dark matter (lower graph, green dashed lines). The monochromatic gamma-ray flux is too 
small(large) in the higgsino(wino) case for explaining the Fermi excess at ~ 130 GeV when 
assuming an Einasto profile. We also show the flux obtained for a boost factor BF = 5.6 
in the higgsino case and when assuming that the wino accounts only for 45% of the dark 
matter density, respectively (red solid lines). In both cases, the monochromatic contribution 
to the flux would yield a good fit to the Fermi data (blue dotted lines). However, the contin- 
uum contribution to the gamma-ray flux overshoots the measured flux at lower energies, and 



therefore both scenarios can be ruled out. The Fermi data are taken from 12 
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3 Constraints from antiprotons 



Constraints that are complementary to the continuum gamma-ray spectrum arise from 
the production of antiprotons in the decay and fragmentation of the weak gauge bosons, 
Higgs bosons or quarks produced in dark matter annihilation or decay The resulting 
primary flux of antiprotons can be potentially observed in the cosmic ray spectrum 
measured at the Earth. Antiprotons are produced with a rate per unit of kinetic 
energy and unit volume at a position r with respect to the center of the Galaxy given 
by 

Qp(E,t) = { n m ™ T xVBR,- 1 for ' (3.1) 



DM 



2w^~ aVt °t I J dE 1 annihilation 



where dN^/dE is the number of antiprotons per kinetic energy and per decay produced 
in the annihilation/decay channels / with branching ratio given by BRy, and p,i m (r) is 
the dark matter density profile. 

The antiproton-to-proton flux ratio observed by the PAMELA satellite experiment 
[36] agrees well with the expectation for the secondary antiproton spectrum produced 
by spallation of ordinary cosmic rays on the interstellar medium. Therefore, upper 
limits on a possible primary component from dark matter annihilation or decay can be 



derived (see e.g. 37-42]). The main source of uncertainty enters due to the propagation 
of antiprotons in the Galactic magnetic field from the production point to the Earth. 
In order to obtain lower limits on the dark matter lifetime and upper limits on the 



annihilation cross section we follow here the approach described in detail in Refs. 40 



42 . In particular, we describe the propagation by a cylindrical two-zone diffusion 



model using three representative sets of parameters leading to minimum, medium and 



maximum antiproton fluxes consistent with the charged cosmic ray abundances 43 . In 
order to obtain conservative limits we use the minimal secondary antiproton spectrum 



from Refs. 40,44 , and use a lower cut of T > 1.5 GeV for the kinetic energy. 

For the case of dark matter decaying into Zu, his, W£ with the same branching ratios 
used above, the 95%C.L. lower limits on the dark matter lifetime is also shown in Fig. [2] 
as blue shaded regions, for the three sets of propagation parameters and the Einasto 
dark matter profile. For a branching ratio much smaller than one, the 95%C.L. lower 
limit on the lifetime is given by r > (0.29, 1.34, 3.31) x 10 27 s for the minimum, medium 
and maximum propagation models, respectively. Compared to the uncertainty arising 
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from the propagation model, the dependence of the lower limit on the dark matter 
profile is rather weak. As can be seen in Fig. [2j the constraints from antiprotons require 
the branching ratio to lie above ~ (9(1%, 5%, 10%), depending on the propagation 
model, if we require that the partial lifetime into monochromatic photons accounts 
for the gamma feature at 130GeV. Even for the minimal propagation model, these 
lower limits on the branching ratio are more stringent than the ones obtained from the 
continuum gamma-ray flux. However, they are also more dependent on astrophysical 
uncertainties related to propagation as well as the secondary antiproton flux Ell. For 
comparison, we also show the lower limits on the lifetime obtained from just requiring 
that the primary antiproton flux does not overshoot the observed antiproton-to-proton 
ratio. The resulting limits are somewhat weaker, r > (0.64, 2.81, 6.62) x 10 26 s. For the 
most conservative case, i.e. for minimal propagation, the resulting lower limit on the 
branching ratio is comparable to the one obtained from the continuum gamma-ray flux 
(see Fig.g. 

For the case of annihilating dark matter into W + W~ , the antiproton constraints 
shown in Fig. [4] correspond to upper limits on the annihilation cross section ov < 
(6.1, 0.69, 0.25) x 10 _25 cm 3 /s when taking the secondary flux into account. The resulting 
lower limits on the branching ratio ov^j ov are more stringent or comparable to the ones 
from the continuum flux for maximum or medium propagation, respectively, as can be 
seen in Fig|4j In the conservative case, when taking only the primary flux into account, 
one obtains rather weak limits ov < (2.8,0.31,0.12) x 10~ 24 cm 3 /s that are comparable 
to the continuum flux only in the maximum case. The antiproton constraints are very 
similar for the case when dark matter annihilates into / yy/ , yZ/WW/ZZ final states (see 
Fig. [5J, except for the antiprotons from 7Z that contribute for very large branching 
ratios. 

Altogether, we find that antiproton constraints are typically the leading ones for de- 
caying dark matter, while the continuum gamma-ray flux is more or equally important 
for annihilating dark matter. 

4 Morphology of the Fermi excess at 130 GeV 

Apart from constraints arising from the gamma-ray spectrum or from antiprotons, 
the Fermi-LAT data also provide valuable information about the spatial distribution 
of the gamma-ray flux observed towards the Galactic center (the angular resolution 
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of the LAT is A9 ~ 0.2° [12]). The hints for a feature in the spectrum are based 
on about 69(57) photons collected within 43 months inside search regions 3(4) (solid 



angle Att/4n = 0.025(0.015)) defined in [12] and the SOURCE sample, and 53(46) 
for the ULTRACLEAN sample. Consequently, any information about the detailed 
spatial distribution of the excess within these regions is at present necessarily limited 
by statistics. Nevertheless, the amount of data allows to test the compatibility with 
certain spatial distributions. In particular, in the analysis of [19] it was found that 
using a Gaussian template the excess appears to originate from a rather small region 
~ 3°FWHM, possibly with a small offset from the Galactic center of 1.5° along the 
Galactic plane. Additionally, it was argued that an Einasto or contracted NFW profile 
yields an even better description, when also taking an offset into account. 

In the following, we leave aside the question about the origin of a possible offset 
in connection with a dark matter explanation. Instead, under the tentative hypothesis 
that dark matter is indeed responsible for the Fermi feature, we discuss the impact 
of the apparent morphology on decaying versus annihilating dark matter. The flux 
of photons arriving at the Earth per unit of energy and solid angle from dark matter 
decay or annihilation is given by the line-of-sight integral over the dark matter density 
or the density squared, respectively, along the direction of observation (see e.g. [46]), 

/ dsp dm (r) decay 

T^m D M i.o.s. 

5 (4.1) 

— 2~^" / - Pdm( r ) 2 annihilation 

m DM l.o.s. ^ 



where r = \/(ro — scos£) 2 + (ssin£) 2 , r$ ~ 8.5kpc is the distance of the Earth from 
the Galactic center, and £ the angle measured with respect to the direction towards the 
center of the Galaxy. As a reference, we consider the density profiles given, respectively, 



by a generalized Navarro-Frenk- White (NFW) [47] or Einasto profile 48 



Pdm(r) oc - 1 — - exp ( — —{r/r s ) aE ) , (4.2) 

(r/r s ) a (l + r/r s y a \ a E J 

with aE — 0.17, a — 1 and scale radius r s = 20kpc. We also consider two strongly 
peaked distributions with a = (1.3,1.7), that may arise from adiabatic contraction. 
All profiles are normalized such that the local density is po = 0.4GeV/cm 3 , except 
for a = 1.7 where we use po = 0.28GeV/cm 3 corresponding to the maximal allowed 
lcr-value when taking observational constraints from microlensing and measurements 



of the Galactic dynamics into account 49 
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Figure 7: Expected spatial profile of the gamma-ray emission from sections along the Galactic 
plane (A£ = 0.5°, |6| < 5°, upper graph) and perpendicular to the plane (Ab = 0.5°,— 5° < 
t < —2°, lower graph) for decaying and annihilating dark matter and two density profiles 
centered around (£,b) = (— 1.5°,0). For illustration, we also show the count maps indicating 



the shape of the Fermi-LAT excess in the energy range 124.7 — 133.4 GeV taken from Ref. 19 . 
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In Ref . pi it has been argued that the morphology of the Fermi excess at ~ 130GeV 
is consistent with dark matter annihilation assuming the Einasto or NFW profile, pos- 
sibly with a slightly larger than one. In order to obtain a similar emission profile for 
decaying dark matter, a much more peaked dark matter distribution would be necessary. 
In order to illustrate this point, we show in Fig. [7] the expected spatial distribution for 
the contracted NFW profiles. For comparison we also show the corresponding shape for 
annihilating dark matter, and the count maps taken from [19] . Qualitatively, it seems 
that even for the strongly contracted profile with a = 1.7 the expected distribution 
for decaying dark matter is too flat, especially for the profile along the Galactic plane. 
One may speculate that, in the perpendicular direction, the excess could be less peaked. 
However, a definite statement is presently impossible due to statistical limitations. 

For decaying dark matter, one would expect also a comparably large signal arising 
from the Galactic halo, where no excess has been observed. Therefore, it is instructive 
to compare the halo limit with the required lifetime to explain the excess in the central 
region. The partial dark matter lifetimes r lv that are required in order to explain the 
feature in the Fermi gamma-ray spectrum present within various search regions defined 
in 12 , and for the various assumptions of the radial dark matter distribution, are shown 



in Fig. [8} Although the best-fit values for the considered search regions are compatible 
with each other, they are in tension with the lower limits on the lifetime obtained by 
the recent analysis of the Fermi collaboration |9|, also shown in Fig. [8] Note that the 
Fermi analysis is based on a rather wide search region, and therefore sensitive to the 
dark matter flux originating in the halo of the Milky Way. Due to the linear scaling of 
the flux with the dark matter density in the case of decaying dark matter, one typically 
expects a larger flux from the halo compared to the center than for annihilating dark 
matter. Thus, for the case of decaying dark matter the gamma-ray line should have 
been observed also within the search region used by the Fermi collaboration, provided 
that the radial distribution of dark matter can be well described by the considered NFW 
or Einasto profiles. The discrepancy is less pronounced for the contracted profiles which 
are more cuspy towards the Galactic center. On the other hand, for annihilating dark 
matter, the cross section required to explain the Fermi gamma-ray feature is compatible 



with earlier upper limits assuming NFW or Einasto profiles 12 



Clearly, this tension disfavors the explanation of the excess in terms of decaying 
dark matter. However, because of the large uncertainty of the actual dark matter 
distribution close to the Galactic center as compared to the halo [2 16 50 , and since 
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Figure 8: Best fit values for the partial dark matter lifetime r™, for a decay into monochro- 
matic photons, obtained for the search regions 2, 3 and 4 defined in Ref. [12] , for various dark 
matter profiles. The error bars show the 95%C.L. regions. Also shown are the lower limits on 
the partial lifetime from the most recent Fermi LAT analysis (the lower error bars include the 
estimated systematic uncertainty) [9], that are based on a wide search region (Halo Limit). 
The upper and lower figures correspond to the two Fermi data samples used in [12]. 
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one cannot exclude large systematic uncertainties at present [9), we believe that it 
would be premature to rule out decaying dark matter as an explanation of the Fermi 
excess based on this tension. 

5 Conclusion 

In light of tentative evidence for a gamma-ray line observed by Fermi-LAT near the 
Galactic center, we have investigated the consequences for several prototype scenarios 
of decaying and annihilating dark matter, motivated by supersymmetric models with 
gravitino, higgsino or wino-like LSP. We find that, independently of the actual dark 
matter distribution, the consistency of continuum and monochromatic contributions 
to the photon spectrum from dark matter decay or annihilation requires a branching 
ratio into monochromatic photons larger than BR 7 > 0.5%. Both higgsino and wino 
dark matter can be ruled out because of a too large continuum flux, independently of 
the production mechanism, while gravitino dark matter with wino NLSP is compatible. 
We have also investigated constraints arising from the primary antiproton flux, which 
depend on the adopted propagation model. The resulting limits are comparable to 
the ones from continuum gamma-rays for annihilating dark matter, but can be more 
important for decaying dark matter. The morphology of the tentative excess, if con- 
firmed, could yield valuable information on the distribution of dark matter close to the 
Galactic center. When assuming conventional NFW or Einasto profile functions, the 
required partial lifetimes for decaying dark matter r 7J , ~ (1 — 3) x 10 28 s are in tension 
with lower limits obtained from the Galactic halo. This tension is less pronounced for 
more cuspy contracted profiles. If decaying gravitino dark matter is indeed responsible 
for the excess, the dark matter density should be enhanced in the Galactic center region 
compared to conventional models. 
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Appendix 



In the following we collect the formulae needed to evaluate the various branching ratios 
of gravitino LSP decays in supersymmetric models with bilinear R-parity breaking. 
The couplings of the gravitino field ijj v to matter and gauge fields fields are given 



by 51 52 



C 



(X7YM^ + c.c.; 



1 



(5.1) 



\/2M v/w ' v pr ' r " ' 4M' 
Here x stands for the left-handed lepton and higgsino doublets /, h u , hd, and for the 
corresponding scalar leptons and Higgs fields /, H u , H^, F vp denotes the fields strengths 
of U (l)y and SU (2) gauge interactions, B up and W* , and A the corresponding gauginos 
b and w 1 . 

In the case of bilinear R-parity breaking the transition from weak to mass eigen- 



states has been worked out in detail in 53 , 54 . Inserting these field transformations 



in Eq. (5.1) and shifting the Higgs fields around their vacuum expectation values, one 



finds for the R-parity breaking gravitino couplings 



Ad 



-{ {^hid^h + iK zi m z Z p ) ViYl^v + iV2K wi m w W eff^^v 



+ im z (Ud P Z u + £^4,) va V^ A 
+ iV^mw^id^W-e^cr^tPx} + h.c. 



(5.2) 



The couplings k and £ are determined by the R-parity breaking parameters Q, the 



gaugino masses M 12 and the \i parameter |53| 54 

c 



K hi 



1 - sin 2(3 
1 + sin 2/3 



m% (M 1 cos 2 6 W + M 2 sin 2 9 W ) 
M x M 2 \i 

m% [M x cos 2 9 W + M 2 sin 2 9 W ) 



MiM 2 /i 



£71 



V2 

A 

4= I 1 -sin2/3 „ , , 

Ci a . a mz (M 2 - Mi) 
-^cos^sm^ , 

Q m z (M 2 sin 2 6 W + Mx cos 2 6 W ) 



y/2 

\ftM 2 



M X M 2 



(5.3) 

(5.4) 

(5.5) 
(5.6) 
(5.7) 
(5.8) 
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where 6w is the weak angle. For a gravitino LSP one has Mi, M 2 , // > ^3/2- Note that 
the gauge boson couplings £ satisfy the relation 



izi + tan ^w^i — iwi 



(5.9) 



The gravitino couplings (5.2) have previously been obtained in an effective Lagrangian 



approach [39 The couplings £ = 0(mz /7713/2) correspond to dimension-6 operators 
and the leading terms of the couplings k to a dimension-5 operator, which gives 



r(^ 3/ 2 -> hui) ~ r(^ 3/ 2 -> ^) ^ 2 r ^ 3 / 2 H/+e * : 



(5.10) 



for mz/ 1713/2 <C 1. The sub leading contributions (^(ml/m 2 ^) to k correspond to 
dimension-7 operators which were not considered in |39|. 



Based on the calculations of 55 , 56 one can now obtain the partial two-body grav- 



itino decay widths by matching coefficients. The result reads 

,3 



r(^3/2 -> hUi) 

T(^3/2 -> IVi) 
r(^3/2 ->■ ZVi) 

r(^ 3/ 2 w+e-) 



m 



384ttM 



3 

3/2 



3/2 2 o4 
2 ^/liP/i > 



TO 



647^M2 |S7^I ' 



3/2 



m 



3/2 



1927rM 2 ' 



J 3/2 

$z ( K wiHw + 16 — — Re («;,„£„, j)GV + 6|^i| 2 -^W 
^3/2 



The functions /3 a , H a , G a and F a (a = h, Z, W) are given by 55 

Ml 



(3a 

G a 
F n 



m 



10 



2 ' 
3/2 

Ml 



+ 



Mi 



2 4' 
m 3/2 m 3/2 



1 Ml 



2m l/2 ' 

2 Af 2 1 Mf 
1 + -^ + 



3 m3 /2 3 m3 /2 



(5.11) 
(5.12) 
(5.13) 

)• 

(5.14) 

(5.15) 
(5.16) 
(5.17) 
(5.18) 



'Note the change in notation compared to 39 . 
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Gravitino Dark Matter 




Figure 9: Branching ratios of two-body gravitino decays for two representative examples. 
Wino NLSP: M 2 = 1.1 rn 3 /2, Mi = \i = 10 m 3 / 2 , and higgsino NLSP: fj, = 1.1 m 3 / 2) 
Mi = 10 m 3 / 2 , M 2 = 1.9 Mi. 

The decay width r(?/>3/ 2 — >■ 7^) is of particular interest since it determines the 
strength of the gamma line at the end of the continuous spectrum. Neglecting correc- 
tions 0(m 2 z /my 2 ), one finds for the branching ratio 



BR^ 3 ;° {M l~y\ (5.19) 
7 2V2G F M\Ml 



where a is the electromagnetic fine-structure constant and Gp is Fermi's constant. 

In Figure [9] the different branching ratios are shown as functions of the gravitino 
mass for the two representative cases of wino and higgsino NLSPs. 
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